Color superconductivity and nondecoupling phenomena in (2+l)-dimensional QCD 
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The possibility of generating a color superconducting state in (2+l)-dimensional QCD is analyzed. 
The gap equation in the leading, hard-dense loop improved, one-gluon exchange approximation is 
£h ' derived and solved. The magnitude of the order parameter is proportional to a power of the coupling 

constant. For an asymptotically large chemical potential, a qualitatively new [with respect to the 
(3 + l)-dimensional case] phenomenon of nondecoupling of the fermion pairing dynamics from the 
infrared one is revealed and discussed. 
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The study of color superconductivity has attracted a great deal of attention in the past few years. The first 
theoretical studies of the subject appeared over twenty years ago |^^]. For many reasons, however, its importance 
was not fully appreciated until recently ||[||. The last couple of years were marked by intensive theoretical studies, 
producing further understanding of the subtleties of physics behind color superconductivity ]5|-|27|] (see Refs. |28 29 
O ,' for a review). 

The color superconducting phase is expected to appear in sufficiently dense quark matter. It is quite likely, therefore, 
that it exists inside some neutron stars ]30| ]. Of course, the hope is that the superconducting order parameter is large 
enough to result in clear observational signatures. While this still remains to be confirmed, more theoretical work is 
• • ■ needed before any solid quantitative predictions would be available. 

In this paper, we study the color superconducting phase in (2 + l)-dimensional QCD at zero temperature and finite 
density. This theory is superrenormalizable and, therefore, asymptotically free. The conventional wisdom is that at 
' asymptotically high quark density the dynamics in asymptotically free theories is weakly interacting and therefore 
drastically simplified. This is the case in (3 + l)-dimensional QCD ||-|^]. One of our goals is to check the validity of 
this picture in (2 + 1) dimensions. We study the problem of color superconductivity using the conventional method 
of the Schwinger-Dyson (SD) equation in the hard-dense loop (HDL) improved one-gluon exchange approximation. 
By solving the SD equation analytically, we obtain an approximate solution for the color superconducting gap. The 
result is proportional to a power of the coupling constant. This contrasts the situation in (3 + l)-dimcnsional QCD, 
where the expression for the gap is nonanalytic in the coupling constant [A ~ fiexp(— C/y/a^) with fi being the quark 
chemical potential and a s = a s (n)] when the coupling a s is weak |^-^|. 

Is the conventional wisdom about weakly interacting dynamics at high quark density in (2 4- l)-dimcnsional QCD 
justified? A detailed analysis of the solution of the SD equation in this theory reveals a qualitatively new [with 
respect to the (3 + l)-dimensional case] phenomenon: nondecoupling of the dynamics of the fermion pairing from the 
nonperturbative infrared one, even in the case of asymptotically large values of the chemical potential. We argue that 
this phenomenon is common in lower dimensional models. The relevance of the 1/iV/ expansion, with Nf the number 
of fermion flavors, for solving difficulties induced by this phenomenon is pointed out. 
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This paper is organized as follows. In Sec. |J, we discuss different types of the Dirac and Majorana masses, paying 
special attention to their transformation properties under parity which itself could be defined in several ways. We 
argue that, in the QCD model at hand, only that color breaking Majorana mass is generated which preserves the 



global SU(2) "pre-flavor" symmetry. The hard-dense loop gluon propagator is derived in Sec. III. Then, in Sec. |IV| , 
we derive the SD equation in the HDL improved rainbow approximation. We solve the equation analytically and 
present an approximate solution. In Sec. |V] the nondccoupling phenomenon taking place for large values of the 
chemical potential in (2 + l)-dimensional QCD is discussed. Our conclusions are given in Sec. VI. In Appendix [A|, 
we describe our notation and list some useful formulas, used throughout the main text. And in Appendix |BJ, we give 
an approximate analytical solution of the gap equation. 

II. ORDINARY AND MAJORANA MASS TERMS 

The (2 + l)-dimensional models are somewhat different from their analogues in (3 + 1) dimensions. In particular, 
the global "flavor" symmetries have somewhat different status. From studying color superconductivity in different 
models, we know that often color breaking is accompanied by locking with flavor symmetries |l2] , |l6| ]. Before studying 
the SD equation for a possibility of color superconducting gap, we discussed some general symmetry properties of all 
possible Dirac and Majorana masses in (2 + l)-dimensional QCD. 

In (2 + 1) dimensions, the irreducible spinor representation is two dimensional. For many purposes, however, it is 
more convenient to work with the four component spinors that contain two different irreducible representations. To 
make a connection with the usual terminology in (3 + 1) dimensions, we could assign a notion of flavor to each four 
component spinor, while using "pre-flavor" for two component irreducible spinors. 

In the case of four component spinors, the kinetic term of quarks, 

£kin = 4> (i^Dn + ip, (1) 
has U(l)v x SU(2) (pre-flavor) symmetry pl|-p4[, generated by 

X, T 1 = 7 3 7 5 , T 2 = rf, T 3 = 7 5 . (2) 
One could also add a Dirac mass term that respects this symmetry, 

toi^TV- (3) 

This mass, however, is odd under parity 

i/j(xi,x 2 ) -> r r 1 4>{-xi,x 2 ), (4) 

and, furthermore, there are no modification of the parity transformation (there are at least three other possible 
definitions, ip — > ^y°j 2 ^p' , ip —* 7 5 7 1 V'' and ip — > 7 3 7 1 ?/>', respectively) that would make the mass term in Eq. (|^) parity 
even. Also, having the mass term as in Eq. (|j|), would lead to the generation of the Chern-Simons term at one loop 
order. 

There are no other mass terms which would respect the U (l)y X SU(2) (pre-flavor) symmetry ^6[. For example, 
each of the following mass terms 

rriQipip, irri2'tpT 2 ip, im^ipT 3 ^, (5) 

breaks the symmetry down to U(l)v X U(l), generated by (X, T 1 ), (X,T 3 ), and (1,T 2 ), respectively. Note that the 
last two terms are even under parity defined in Eq. (^). Although no definition of parity exists that would be respected 
by all three mass terms in Eq. (||), a different definition of parity (ip 7 5 7 1- 0' or ip — > 7 3 7 1 ^'') could make the first 
term in Eq. (Q) parity even. Also, since the SU (2) transformations could change each term in Eq. (j^) into another, 
all three of them are physically equivalent. 

As a result of spontaneous breaking of the color symmetry, different kinds of Majorana masses could be generated. 
In analogy with the ordinary masses, the following four structures are possible: 

Ai/S c V, U(l) by T 2 , (6a) 

AxV^TV, U(l) by T 3 , (6b) 

A 2 ^ c T 2 iP, SU(2) by r\T 2 ,T 3 (6c) 

A 3 ^ C T 3 V, 17(1) by T 1 , (6d) 
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where ip = Cip T denotes the charge conjugate spinor, and C is a charge conjugation matrix, defined by C 1 7 M C = 
— 7J and C — — C T . In Eqs. (|6a| ) - (6d) we also indicated the symmetry of the Majorana mass terms and the 



generators of the corresponding groups. Notice that the original U(l)v is always broken. 
The parity properties of the Dirac and Majorana mass terms are summarized in Table I. 

TABLE I. Parity properties of Dirac and Majorana mass terms. 
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While the original field ip belongs to the fundamental representation of the (pre-flavor) SU(2) group, the charge 
conjugate field ip c belongs, in general, to the anti-fundamental representation. It is clear than that the diquark fields 
such as the Majorana mass terms in Eq. (JsJ) , could belong either to the (antisymmetric) singlet representation (A2), 
or to the (symmetric) triplet one (A, Ai and A3). 

Since the most attractive diquark channel is antisymmetric (e.g., antitriplet in the case of SU(3) color group), it 
becomes obvious that the color symmetry breaking order parameter would be SU(2) singlet as that in Eq. (pq). Note 
that this order parameter is odd under the parity defined in Eq. (^), but even under the parity P3 (see Table I). 



III. HARD DENSE LOOP APPROXIMATION FOR THE GLUON PROPAGATOR 



In dense media, the polarization effects play very important role. Below we derive the polarization tensor in the 
so-called hard dense loop (HDL) approximation. This approximation corresponds to dealing with semi-hard (external) 
gluon lines, a <C k <C /J., and hard (internal) fermion lines. Note that there are no one loop corrections with internal 
gluon and ghost lines in HDL approximation. Our notation and some formulas used in this section are given in 
Appendix 0. 

The analytical expression for the polarization tensor reads 



n 



[iv AB 1 t \ _ 



(k) = AiiraNf 



dqod 2 q 
(2ir) 3 



tr [^T A G{q)YT B G{q + k)} , 



(7) 



where a is a (dimensionful) coupling constant of (2 + l)-dimensional QCD and Nf counts the number of the four- 
component quark flavors, not the pre-flavors. 

In derivation of the HDL polarization tensor, it is sufficient to consider only massless quarks. The corresponding 
propagator is given by 



G(°)(g) = i 7 ° 



f 



50 



is 



A« + - A<-> 

qo - e<? + is sign(e 9 ) 



(8) 



where = \p\ ± /j, and the projectors Ap are defined in Eq. ( |A12| ) in Appendix [a|. Note that e is a small positive 
constant that prescribes the proper integration contour around the poles on the real qo _ axis. At the end of calculation, 
this parameter should be taken to zero. 

By substituting the quark propagator into Eq. (Q), we arrive at 



W vAB {k) ~ -2iiraN f 6 



AB 



dq d 2 q 


tr 


7 M 7 o A (-) 7 , 7 o A (-r 




(2n) 3 


. (qo ~H +i£ sign(e 9 )) (q + k - e q+k 


+ is sign(e 


i+fc)) 



tr 


7 M 7 o A (+) 7 , 7 o A (-r 


tr 


7 M 7 o A (-) 7 , 7 o A (+r 




(go + 4 ~ i£ ) (go + h - e q+k + ie sign(e g+fc )} 


(go -£g + is sign(e,j )) (q + k + e+ +fc - is) 



(9) 
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Here we dropped one term in the integrand that contains its both poles in the upper half of the q$ complex plane. 
Such a term does not give any contribution after the integration over go is performed. 

In the HDL approximation (|fc| <C \q\ ~ fi), the traces appearing in the last expression are simple (see Appendix |a|): 



tr 



tr 



tr 



r7 o A (-y 7 o A M j „ tr [yy A (-y 7 o A (-) 



= 2 



.1/0 



7 , 7 o A ( + ) 7 , 7 o A ( ; ) 



tr 



tr 



7 M 7 o A (+) 7 , 7 o A (- 



-2g» v + 2g^g M - 2 



7 M 7 o A (-) r o A (+)j = _ 2 ,u + 2 M o ,o _ 2 ^ r 
9 9 J MM 



(10) 

(11) 
(12) 



In calculation of H fl ' y (k), we first perform the integration over q using the residue theorem. The result reads (here 
we drop the overall 5 AB ) 



W{k) 



,o , <?\ 0(»-\q\)6(\q + k\-fi) d(W-n)0(p-\q + k\) 



M/ V \q\-\q + k\+ k +ie \q\ - \q + k\ + k - ie 



aNf Ws^-^v°+Sft 



\q + k\ - ,i) 



\Q\\g\J \\q\ + \q + k\-ko- is \q\ + \q + fc| + fc - ie 



(13) 



While the first term in this expression is finite, the second term contains a linear ultraviolet divergency. As usual, we 
drop the ultraviolet divergency from the polarization tensor in the HDL approximation (it would go away after the 
the same renormalization procedure as in the case of zero chemical potential) . 

The leftover angular integration could be most easily done by going to a fixed coordinate framework [e.g., such 
that fc — (|fc|,0)]. By doing so, we eventually arrive at the following result for the polarization tensor in the HDL 
approximation: 



n 00 (fc ,fc) = n,(fco,fc), 

n 0i (k a ,k) = k ^-u l (k ,k), 
\k\ 2 

U ij (k ,k)= 



k l W 

W 



n t (fc ,fc) 



^J_n ; (fc ,fc), 



(14) 
(15) 

(16) 



where 



7T/2 



n ; (fco,fc) = ^^ / 



\k\ cosip 



\k\ cos if 



= M 2 e{k 2 )\^-l-i9{-k 2 



n 

tt/2 



n t (fc ,fc) = ^^ / d(p 



fco — |fc| cos ip + ie kg + \k\ cos <p — ie 

' l g l TO ^ sin2 ^ IfclcosysinV +Wy | =M 2__ n; ,, || 

fc — \k\ cos ip + ie k + \k\ cos <p — ie J |fc| 2 



A: 2 



fc 2 



(17) 



(18) 



with fc 2 = fco — |fc| 2 . Here we use the notation M 2 = 2afiNf for the Debye mass. We should point out that the 
expression for H;(fco, fc), as written in Eq. (p"7j), is valid only for real values of fco. The expression for n^fco, fc) as a 
function of complex fco has the following form: 



n ; (fco,fc) = M 2 



1 



l-(|fc|/fc ) 2 



- 1 



(19) 



Our expressions in Eqs. ( |17|) and (|18|) resemble the results in the (3+ l)-dimensional case ||37j— |3S[] . The corresponding 
polarization tensor has a nonzero imaginary part for space-like gluon momenta. This is related to the so-called Landau 
damping of the gluon field. 

The gluon polarization tensor is transverse 
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k^(k Q ,k)=0, (20) 

and the following representation holds: 

n M „ = -o$n, + o$ (n t - a/ 2 ) . (21) 

Here we make use of the projectors of the magnetic, electric and longitudinal gluon modes introduced in Ref. 0, 

0« =g^- u »u v + ^, (22a) 
|fc| 2 

00=^-^-^, (22b) 

0( 3 ) = ^, (22c) 

with u M = (1, 0, 0, 0) and = k^ — (u ■ k)u^. 
Then, we rewrite the inverse gluon propagator as 

(v AB (k , k)) 1 = iS AB (fc 2 - n t ) o« + i5 AB (fc 2 + n t - m 2 ) o$ + l s AB ^o^, (23) 

where £ is a gauge parameter. Then, by making use of the properties of the projection operators, we invert this 
expression and arrive at the final form of the gluon propagator in HDL approximation 

K B (^ k) = -#^o« - ^ fc2 + n t _ M2 off - w^^o®. (24) 

By using the conventional assumption that the pairing dynamics is dominated by the momenta |fc | <C |fc|, we arrive 
at the following approximate expression of the propagator in Euclidean space (feo = ik^): 

W^iiki, fc) ~ — S AB — ^ 0$ - S AB i 0% - 5 AB ^O^. (25) 

^ y ' |fc|3 + M 2 |fc 4 | M kj + \k\ 2 + M 2 kj + \k\ 2 ^ V ' 



IV. SCHWINGER-DYSON EQUATION 



The mechanism of breaking of color symmetry in dense QCD is well known. The dynamics of quasiparticles around 
the Fermi surface is affected by the famous Cooper instability. This causes the perturbative vacuum to rearrange 
so that an energy gap is formed in the spectrum of quasiparticles. This also is accompanied by the appearance of 
a color antitriplet (in the case of three colors) condensate that "breaks" the original gauge symmetry SU (3) c down 
to SU(2) C subgroup. As we discussed in the preceding section, the corresponding order parameter leaves the global 
SU(2) pre-flavor symmetry intact. 

In passing we note that because of the Higgs mechanism, five out of the original eight gluons become massive. 
The value of the corresponding mass functions at zero momentum is expected to be of order afi. This might suggest 
that the mentioned five gluons are rather inefficient in providing interaction in the diquark pairing. If so, one should 
take this into account when studying the SD equation. However, we assume that the Meissner mass function of 
gluons quickly vanishes (approaching the HDL approximation) as we go to the region of momenta larger than the 
superconducting gap. This is what happens in the (3 + l)-dimensional model JtJ, and we assume that this should be 
still valid in (2 + 1) dimensions. 

Our strategy is to use the HDL improved rainbow approximation for the SD equation and then to check whether 
it is reliable [as it happens in the (3 + l)-dimensional case] for an asymptotically large chemical potential. The SD 
equation written in this approximation is: 

(G(p))- 1 = (g(°»(p)) _1 +4™ J ^rXy ( T q _^ A) t ) <W ( ^ _ (T °4 )T ) 2V(?-P), (26) 
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where V tiv {k) is the propagator of gluons with the overall factor S AB omitted. We denote the full quark propagator 
by G(p), and the perturbative one by G^°\p). 

As is argued in Sec. ||, the order parameter should be a color anti-triplet and a pre-flavor singlet, e s,ah {^T 2 ^^ . 
Therefore, the inverse of the full quark propagator should have the following general structure (neglecting the wave 
function renormalization): 



(Gwr 1 = -i ( 



A, 



(27) 



where A p = 7 °At 7 ° and (A p ) afc = e ab ^ (A; A; + A+A+). 

By inverting the expression in Eq. (]27j), we arrive at the following propagator: 



G(p) 



Ri S 

S i?2 



where 
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R 2 = l a ° 
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^_ (e -)2_| A -|2 p 2_ (£ + )2 _| A +| 2 



-x 27 ° 
-x 27 ° 



A; 



A 



A; 



PO + £p 

A: 



Po 



Po 



A-A+ 



A+A; 



p2_ (e -)2_| A -|2 p 2_ (£ + )2 _| A +p 



(28) 

(29) 
(30) 
(31) 



and (2i) (, = S ab - dffi, (I 2 )ab = S^. 
Thus, the gap equation reads 

.8ira f dq 4 d 2 q 



A p A p 



A + A + = -i- 



(27T) 



rA-A+r 



^A+A-r 



g 4 2 + (6 9 -)2 + |A 9 -p g| + ( e +)2 + |A+|2 



T>^(iq 4 - ip 4 ,q- p). (32) 



Among two gap functions, A~ and A+, only the first one plays the role of the gap in the quark spectrum around the 
Fermi surface. In addition, the main contribution in the right hand side of Eq. (^) comes from the first term also 
proportional to A~. As a result, the equation for A~ approximately decouples, 



Ana f dq 4 d 2 q A"tr [yA+^A"] 



£V(% - ipi,q-p)- 



(33) 



(27T)3 g| + ( e -)2 + | A -|2 »» 

By taking the trace and performing the angular integration (see Appendix |a| for details), we arrive at the following 
approximate equation: 



a 

()TT 

+ 



dq 4 dqA Q 



g | + (e-)2 + | A -|2 
2 



^(q-p) 2 + (M 2 \q 4 ~p 4 \) 2 ^ 



^{q-p) 2 + {q 4 -p 4 ) 2 +M 2 V(g-p) 2 + fe-P4) 2 



(34) 



To get a analytical estimate for the solution, we assume that the gap is only a function of the time component of 
the momentum, i.e., A~(p 4 ,p) ~ A~(p 4 ,/i). Then, we can easily perform the integration over e~ = q — fi. The 
approximate result is 



A"(p 4 )^ ~ 



d(? 4 A (g 4 ) 



|A 9 "| 2 



1 



(Af 2 |(74-p 4 |) 1 / 3 M \q 4 -p 4 \ 



(35) 



Because of the approximations, this last equation makes sense only with a finite ultraviolet cutoff which should be 
chosen at the scale of M. 

Now by making use of the same method as in Ref. [Q and using the Landau gauge (£ = 0), we straightforwardly 
obtain an analytical solution (see Appendix p] for details): 



G 




A -(p 4 ) = c|A -| J t , , c «i. 6j (36) 



with the overall normalization constant c chosen so that A (P4)\ P4 ~^- \ — |A |. The value of the gap itself is equal 
to 

i a °i*4 = <V (37) 

where b w 0.331. Note that we neglected the effect of electric gluon modes because their contribution is strongly 
suppressed compared to the effect of the magnetic modes. 

How reliable is the HDL improved rainbow approximation in (2 + l)-dimensional QCD? Because of a singular 
nature of the gauge fixing term in Eq. j35|), we see that the above result for the gap is subject to a large gauge 
dependent correction in a general covariant gauge. This indicates that the HDL improved rainbow approximation 
may not be reliable. In the next section, we will show that the origin of this problem in (2 + l)-dimcnsional QCD is the 
phenomenon of nondecoupling of the dynamics of the fermion pairing from the nonperturbative infrared dynamics. 
This phenomenon is intrinsic for lower dimensional field theories and does not take place in dense QCD in (3 + 1) 
dimensions. We will also discuss a way of solving this problem. 



V. NONDECOUPLING PHENOMENA IN THE GAP DYNAMICS IN (2 + l)-DIMENSIONAL QCD 

Let us first recall that the HDL improved rainbow approximation is reliable in (3 + l)-dimensional dense QCD at 
asymptotically high quark densities. This is related to the fact that a large chemical potential /x provides a large 
value of the gap in (3 + 1 ) dimensions, namely, | Aq | ^> Aq C d for a sufficiently large /i |^-^ . This is enough to justify 
the use of the one gluon exchange approximation for the kernel of the SD equation. Indeed, the integration over the 
gluon momentum in the SD equation is dynamically cut at a momentum k ~ |Aq | 3> A-qcd in infrared, insuring that 
the dominant contribution is given by the hard gluons only. As a result, the diquark pairing dynamics is completely 
decoupled from the strong infrared dynamics. 

The situation in dense QCD in (2 + 1) dimensions is different. While the contribution of hard fermion loops still 
dominates in the polarization operator for large /u,, the value of the gap (|37| ) obtained in the HDL improved rainbow 
approximation in the Landau gauge is small. Therefore, a large region of soft gluons contributes to the SD equation 
(moreover, it dominates). It is clear therefore that the one gluon exchange approximation is not reliable in this case. 

What is the physical origin of such a drastic difference between the (3 + 1)- and (2 + l)-dimensional cases? We 
believe it is a diminishing role of the Fermi surface in the fermion pairing dynamics in lower dimensions. Indeed, the 
main effect of the Fermi surface in the fermion pairing is reducing the initial dimension d + 1 to the effective dimension 
1 + 1. This is crucial for d = 3. However, for example, for d = 1 there is no such a reduction at all: there is 1 + 1 
dimensional dynamics from the outset The present case with d = 2 is intermediate, and our analysis shows how 
different it is from the d — 3 case. 

The importance of the region of soft gluon momenta in the SD equation implies that the pairing dynamics in 
(2 + l)-dimensional QCD is intimately connected with the strong infrared dynamics. In fact, there is no decoupling 
of these two dynamics even for asymptotically large values of the chemical potential. 

Let us show that this point is directly related to the existence of an upper bound for the value of the gap in the 
(2 + l)-dimensional QCD. In particular, |A^~| cannot be much larger than a. Indeed, if the value of the gap were 
much larger than a, then the HDL improved rainbow approximation would be reliable. But the only solution of 
the corresponding equation is given in Eq. (j37|), showing that |Aq | <C a. Thus, the assumption of large |Aq | is 
self-contradictory. 

Therefore in order to solve the problem for large /i in (2 + 1) dimension, one should put the non-perturbative infrared 
dynamics under control. Does it imply that expression (|37]) obtained in the HDL improved rainbow approximation 
is irrelevant? We do not think so: we will argue that it can be qualitatively correct in the case of a large number of 
fermion flavors Nf. 



1 In passing, we note that in the 1 + 1 case, in the HDL approximation, the Debye screening mass is proportional to a?2 (where 
«2 is the gauge coupling) and is independent of /j, at all! The Debye mass actually coincides with the famous Schwinger mass, 
appearing for fi—0 in 1 + 1 [see also a discussion in Ref. jicft]. This observation certainly implies the breakdown of the HDL 
expansion in 1 + 1 dimensional QCD. 
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The fact that the 1/Nf expansion can be useful in (2 + l)-dimensional QED and QCD was recognized long ago 
[fj3| j36f . The point is that because these models are renormalizable theories, they are asymptotically free for any 
number Nf. Though being very nontrivial, the 1/iV/ expansion is helpful in putting under control the infrared 
dynamics. The crucial point is the selection of a "right" gauge in the leading order in 1/Nf. In particular, appropriate 
Ward identities have to be satisfied in that gauge. In other gauges, the result can be found by gauge-transforming 
Green functions from the "right" gauge to those gauges. The description of this very nontrivial procedure can be 
found in Ref. gl). For our purposes, it is enough to know that though the "right" gauge is different from the Landau 
one, the results obtained in the improved rainbow approximation in the latter are qualitatively reliable: its special 
role amongst other covariant gauges is connected with the approximate validity of the Ward identities in the improved 
rainbow approximation using the bare vertices. 

It is reasonable to assume that the Landau gauge is special in our problem either. Indeed, let us recall that, because 
of the Miessner effect, the quark gluon vertices should necessarily receive non-perturbative pole corrections J42j. This 
is a rather general and a model independent property that is related to the Ward identities for the vortex 

function. Following Ref. [42]] , one can show that in a general covariant gauge the contribution of those pole terms 
in the SD equation is large, thus destroying the validity of the approximation with bare vertices. The solution to a 
similar problem was suggested long time ago by Cornwall and Norton p^] . The idea is to choose the gauge where the 
approximation of the SD equation with bare vertices is self-consistent. In their model, it was the Landau gauge. With 
this gauge choice, it was possible to show that potentially dangerous contributions cancel from the SD equation. This 
is possible due to the fact that the pole contribution to the vertex is pure longitudinal. The Landau gauge is special 
because the gluon propagator is completely transverse. Thus, when used in the SD equation, such gluon propagator 
annihilates the dangerous infrared poles. 

Using the same arguments, we see that the Landau gauge is the best gauge among the general covariant gauges 
in the problem at hand too. We should note, however, that our situation differs from that in Ref. 0|. The Lorentz 
symmetry is broken in dense QCD by a nonzero chemical potential. As a result, the pole structure of the quark gluon 
vertex is proportional to p M — g** = (p° — q°, (p — g)/3), rather than p^ 1 — q 1 * as in Ref. f44| . While the four- vector 

— q^ is indeed annihilated by the magnetic modes of gluons, it is only partially annihilated by the electric ones. It 
is fortunate that the electric gluon modes are screened much stronger and their effect on diquark pairing is negligible 
in the model at hand. So, we could still justify the use of the Landau gauge, in which all most dangerous infrared 
contributions are canceled from the SD equation. 



These arguments lead us to believe that expression (37) can be qualitatively reliable for large Nf. We rewrite the 
result in the following form: 

(38> 

where a = Nfa and the constant b is of order one. We got this estimate by making use of conventional large Nf 
rescaling of the gauge coupling. Though a color-flavor structure of the condensate can vary with TVj, we expect that 
the magnitude of the gap should be similar for different, but same order, large Nf. 

VI. CONCLUSION 

We studied the dynamical generation of the color superconducting order parameter in (2 + l)-dimcnsional QCD at 
asymptotically large baryon density and zero temperature. For this purpose, we used the conventional method of the 
Schwinger-Dyson equation in the HDL improved rainbow approximation. It was found that the order parameter in 
this theory is similar to that in the S2C phase of the (3 + l)-dimensional two flavor QCD. 

By solving the SD equation analytically, we obtained an approximate solution for the color superconducting gap. 
The parametric dependence of the result on the coupling constant is given by a power law. This is qualitatively 
different from the situation in the (3 + l)-dimensional QCD, where the expression for the gap is exponentially small 
for small values of the coupling constant . 

The analysis of the solution of the SD equation revealed a qualitatively new [with respect to the (3 + l)-dimcnsional 
case] phenomenon: nondecoupling of the dynamics of the fermion pairing from the nonperturbative infrared one, even 
in the case of asymptotically large values of the chemical potential. We believe that this phenomenon is common in 
lower dimensional models and deserves further study. 

It would be also interesting to study the interplay between the dynamics of color and chiral condensates in (2 + 1)- 
dimensional QCD. It is known that, at zero quark density a chiral condensate can occur only if the number of quark 



flavors is less than a critical value Ny , which is a function of the number of colors N c p5 36 1. An interesting question 
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is whether, in this case, there is a phase transition which restores chiral symmetry at a critical value of the chemical 
potential and how the character of this phase transition depends on the values of Nf and iV c . 

In the analysis of the SD equation we assumed that the Meissner effect is mostly irrelevant for the pairing dynamics 
in dense quark matter. While this was well justified in (3 + l)-dimensional QCD at asymptotic densities j7|fl^5|, 
it remains an open issue in the (2 + l)-dimcnsional case. Indeed, our present analysis suggests that the dominant 
contribution to the gap equation comes from the region of momenta not much larger than |A _ |. The Meissner 
effect, therefore, could play a significant role in modifying the gluon interaction in the diquark channel. It would be 
interesting to study this problem in more detail. 

It is interesting to point out that, because of close similarity of the model studied here with that of the (3 + 1)- 
dimensional two flavor QCD, one should also expect the appearance of five light (Md q <C |Ao|) pseudo-Nambu- 
Goldstone diquarks in the low energy spectrum 42 ] . The argument for the existence of such states is essentially the 
same as in the (3 + l)-dimensional case, see Ref. [ 42] . The only difference should appear in the hierarchy of the scales. 

At last, we would like to mention that some relativistic (2 + I)-dimcnsional gauge models are used by some theorists 
for describing high temperature superconductivity in cuprates Q. It would be interesting to see whether the color 
superconductivity could anyhow be related to the ordinary superconductivity in planar systems. 
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NOTE 



While writing our paper, we learned that a partially overlapping study is being done by Prashanth Jaikumar and 
Ismail Zahed. 



APPENDIX A: NOTATION AND USEFUL FORMULAS 



We use the four-dimensional representation of the Dirac 7-matrices. In addition to the three matrices 7 M (/i = 0, 4, 2) 
there are two more matrices that anticommute with 7^ as well as between themselves, 7 3 and 7 s [here we use the 
same notation for the 7-matrices as one would use in the (3 + 4)-dimensional case]. 

While studying the color superconductivity in the model, it is convenient to introduce the charge conjugate spinors 
as follows: 

^ c (x) = Cj T (x), ${x) = - ty c (x)] T C\ (Al) 
F(x) = -^ T (x)C\ <P(x) = C [4> c (x)f , (A2) 

where C is a unitary charge conjugation matrix, defined by C -1 7 M C = — 7J and C = —C T . 
Note that in momentum space, 

^ c {p) = Cj T (-p), ^p) = -[^(- p )] T c\ (A3) 
j c (p) = -V{-P)C\ HP) = C $ C (-P)] T . (A4) 
The Lagrangian density of QCD could be written as 

£qcd = \$ (jf + M7°) ^ + \i> Cl {i> - M7°) ^ C + \^KV - \$ C -fAZl> C - ^Tr (P^) . (A5) 
The eight component Majorana spinor is defined by 
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1 



( v 



V2 \ 

In this new notation, the bare fermion propagator reads 

i> + /i7° 



(g(°)(p)) _1 



i> - /i7° 

Note the following convenient representations for the diagonal elements of this propagator: 



^ + /i7° = 7° 



(po-e;)A( +) + (p + e+)A(-) 
J* - M7° = 7° \(P0 - e+)A« + (p + e£)A<-> 



(jf + MT )" 1 =7° 



Ai+> 



1 



M7 



1 



-Al+) 



where er 1 = \p\ ± /i and the "on-shell" projectors are 



A(±) = - 
p 2 



1 ± 



a ■ p 



Po - e P 
1 

po - <4 



a = 77. 



-Ai" 



The bare quark- quark- gluon vertex, in the eight component notation, is also a matrix, 



r 







-{T A ) T 



(A6) 

(A7) 

(A8) 
(A9) 

(A10) 
(All) 

(A12) 
(A13) 



While working with the SD equation in Sec. |IV|, we encounter a few Dirac traces. Here are the results for most general 
ones: 



tr 



Yk [ p ) l v k [ q' ) ] = 3^(1 + ee't) - 2eeyV°£ + ee' 



tr 



7 ^7°A^7 l '7 A( e ') = -.9^(1 - ee't) + ( ^° - e'^ ) <T - 



IpI 



(A14) 



(A15) 



where e, e' = ±1, t = cos (p is the cosine of the angle between two- vectors q and p, and irrelevant antisymmetric terms 
are denoted by ellipsis. 

By making use of the traces above, we calculate the following expressions which appear in the SD equation: 

0« tr [VAWyAMl = (1 - t) 2 (Alto) 



q 2 + p 2 — 2qpt ' 
tr ^A^yA^l = (1 + t) - 3 (q ~ P) l^ + (1 - t 2 ) 



2qp 



tr fyA(+yA(-)l =(l + t) 



q 2 +p 2 - 2qpt q 2 + p 2 - 2qpt + (g 4 — p 4 ) 2 ' 

(g-p) 2 + ( g4 -p4) 2 
g 2 + p 2 - 2gpi + (g 4 - Pi) 2 ' 



(A16b) 
(A16c) 



where q = \q\, p = \p\, 94 = — iqo and p^ = — ipo- Making use of these expression, we perform the angular integration 
that appears in the SD equation, 



1 



y/{q-pf + {M 2 Lj) 2 l 3 ^{q - p) 2 + w 2 + M 2 



\/(q-p) 2 +lj 2 



(A17) 



Here M 2 = 2afj,Nf and lu = |q 4 — p 4 \. Although we performed the angular integration exactly, we dropped all 
subleading terms on the right hand of Eq. ( Al7| ), assuming that both q and p remain in the vicinity of the Fermi 
surface. 
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APPENDIX B: APPROXIMATE SOLUTION OF THE GAP EQUATION 



In this Appendix we present an approximate analytical solution to the gap equation (|35|). In order to reduce the 
integral gap equation to a differential one, we introduce an infrared cutoff in the integral at the scale of |Aq |, and we 
approximate the kernel in the integrand by its asymptotes: 



1 



1 



(M2|p 4 |)l/3g 4 • 



yJqj + \Aq\*( M2 ^-P^ 1/3 I M*/3|, 4 |Vs 



for 
for 



We then arrive at the following equation: 



A-(p 4 



dq 4 



I A" | 94 



A" ((74 



,1/3 



dq 4 



94 



4/3 



94 < Pi, 
94 > Pi- 



A-(g 4 ), 



where 



a 

M 5 ' 



(Bl) 



(B2) 



(B3) 



It is straightforward to check that the integral equation is equivalent to the following second order differential equation 



d 2 A-{ Pi ) 4 dA-( Pi ) 



dpj 



j a 3p 4 dp 4 

along with the infrared and ultraviolet boundary conditions: 

dA~(p 4 ) 



I/V3 



A-(p 4 )=0, 



dp 4 

dA~(p 4 ) 

3p 4 — J hA (p 4 ) 

dpi 



P4=|A 



0, 



= 0. 



P4 — OO 



The general solution to the differential equation (B4) is given in terms of Bessel functions, 

V 6 / .1/6 \ / i a — i \ 1/6 



A~(P4 



I An 



P4 



Ji 2 



^6 

1/6 
P4 , 



C 2 



Pi 



Vi 2 



P 4 /6 , 



(B4) 

(B5) 
(B6) 

(B7) 



In order to satisfy the ultraviolet boundary condition in Eq. (B6), we m ust choose C2 = 0. In passing we note, however, 
that when a finite ultraviolet cutoff is used in the gap equation (B2) (one might choose, for example, q 4 — M), the 
integration constant C2 will be nonzero, 



C-2 



2 \mJ 



2/3 



(B8) 



As is easy to check, the effect of this nonzero constant is negligible in the leading order approximation. 
Now, by satisfying the infrared boundary condition, we arrive at the following equation: 



Ji 2 



,1/6 



|A -|i/ 6 



1/6 



J 2 2 



^1/6 

IaTF 



1/6 



Jn 2 



,1/6 



IAqTI 1 / 6 



(B9) 



This equation has infinitely many solutions, which correspond to different solutions to the gap equation. We must 
choose a single solution which gives the gap function without nodes. As one could check, this is also equivalent to 
choosing the solution with the largest value of the gap. Thus, we obtain |Ajj"| = bv where b w 0.331. 



[1] B.C. Barrois, Nucl. Phys. B129, 390 (1977); S.C. Frautschi, in "Hadronic matter at extreme energy density", edited by 
N. Cabibbo and L. Sertorio (Plenum Press, 1980). 



11 



D. Bailin and A. Love, Nucl. Phys. B190, 175 (1981); Nucl. Phys. B205, 119 (1982); Phys. Rep. 107, 325 (1984). 
M. Alford, K. Rajagopal, and F. Wilczek, Phys. Lett. B422, 247 (1998). 
R. Rapp, T. Schafer, E.V. Shuryak, and M. Velkovsky, Phys. Rev. Lett. 81, 53 (1998). 
R.D. Pisarski and D.H. Rischke, Phys. Rev. Lett. 83, 37 (1999). 
D.T. Son, Phys. Rev. D 59, 094019 (1999). 

D.K. Hong, V.A. Miransky, LA. Shovkovy, and L.C.R. Wijewardhana, Phys. Rev. D 61, 056001 (2000). 
T. Schafer and F. Wilczek, Phys. Rev. D 60, 114033 (1999). 
R.D. Pisarski and D.H. Rischke, Phys. Rev. D 61, 051501 (2000). 
S.D.H. Hsu and M. Schwetz, Nucl. Phys. B572, 211 (2000). 

W.E. Brown, J.T. Liu, and H.-C. Ren, Phys. Rev. D 61, 114012 (2000); ibid. 62, 054016 (2000). 

M. Alford, K. Rajagopal, and F. Wilczek, Nucl. Phys. B537, 443 (1999); M. Alford, J. Berges, and K. Rajagopal, Nucl. 
Phys. B558, 219 (1999). 

T. Schafer and F. Wilczek, Phys. Rev. Lett. 82, 3956 (1999); Phys. Rev. D 60, 074014 (1999). 
V.A. Miransky, LA. Shovkovy, and L.C.R. Wijewardhana, Phys. Lett. B 468, 270 (1999). 
LA. Shovkovy and L.C.R. Wijewardhana, Phys. Lett. B 470, 189 (1999). 
T. Schafer, Nucl. Phys. B575, 269 (2000). 

D.H. Rischke, Phys. Rev. D 62, 034007 (2000); ibid. 62, 054017 (2000). 

G. W. Carter and D. Diakonov, Nucl. Phys. B582, 571 (2000); R. Casalbuoni, Z. Duan, and F. Sannino, Phys. Rev. D 62, 
094004 (2000). 

R. Casalbuoni and R. Gatto, R. Casalbuoni and R. Gatto, Phys. Lett. B 464, 111 (1999); D.K. Hong, M. Rho, and 
I. Zahed, Phys. Lett. B 468, 261 (1999). 

D.T. Son and M.A. Stephanov, Phys. Rev. D 61, 074012 (2000). 

M. Rho, A. Wirzba, and I. Zahed, Phys. Lett. B 473, 126 (2000); M. Rho, E. Shuryak, A. Wirzba, and I. Zahed, Nucl. 
Phys. A 676, 273 (2000). 

D.K. Hong, T. Lee, and D.-P. Min, Phys. Lett. B 477, 137 (2000); C. Manuel and M.G.H. Tytgat, ibid. 479, 190 (2000). 
K. Zarembo, Phys. Rev. D 62, 054003 (2000). 

S.R. Beane, P.F. Bedaque, and M.J. Savage, Ph ys. Lett. B 483, 131 (2000). 
D. H. Rischke, D. T. Son and M. A. Stephanov , |hep-ph/0011379| . 
K. Rajagopal and F. Wilczek , |ticp-pli/0012039| . 
V. H. Nguyen and X. Pham, frep-ph/0101332r 
K. Rajago pal and F. Wilcze k, |ticp-pli/0011335 . 
M. Alford, fiep-ph/0102"047| . 

M. Alford, J. A. Bowers and K. Rajagopal, |hep-ph/0009357 . 
R. Jackiw and S. Templeton, Phys. Rev. D 23, 2291 (1981). 
S. Deser, R. Jackiw and S. Templeton, Ann. Phys. (N. Y.) 140, 372 (1982). 
R. D. Pisarski, Phys. Rev. D 29, 2423 (1984). 

T. W. Appelquist, M. Bowick, D. Karabali and L. C. Wijewardhana, Phys. Rev. D 33, 3704 (1986); ibid. D 33, 3774 
(1986). 

T. W. Appelquist, D. Nash and L. C. Wijewardhana, Phys. Rev. Lett. 60, 2575 (1988); D. Nash, ibid. 62, 3024 (1989). 
T. W. Appelquist and D. Nash, Phys. Rev. Lett. 64, 721 (1990). 

H. Vija and M.H. Thoma, Phys. Lett. B 342, 212 (1995). 
U. Heinz, Annals Phys. 168, 148 (1986). 

C. Manuel, Phys. Rev. D 53, 5866 (1996); G. Alexanian and V. P. Nair, Ph ys. Lett. B 390, 3 70 (1997) jhep-th/9609011 |. 
B. Park, M. Rho, A. Wirzba and I. Zahed, Phys. Rev. D 62, 034015 (2000) |hep-ph/9910347| . 

K.-I. Kondo and H. Nakatani, Mod. Phys. Lett. A 5, 407 (1990); C. J. Burden and J. Praschifka, and C. D. Roberts, Phys. 
Rev. D 46, 2695 (1992); P. Maris, ibid. 52, 6087 (1995); P. Maris, ibid. 54, 4049 (1996); V. P. Gusynin, A. H. Hams, and 
M. Reenders, Phys. Rev. D 53, 2227 (1996); I. J. R. Aitchison, N. E. Mavromatos, and D. McNeill, Phys. Lett. B 402, 
154 (1997). 

[42] V.A. Miransky, LA. Shovkovy, and L.C.R. Wijewardhana, |hep-ph/0003327| ; Phys. Rev. D 62, 085025 (2000); ibid. D 63, 
056005 (2001). 

[43] R. Jackiw and K. Johnson, Phys. Rev. D 8, 2386 (1973). 
[44] J.M. Cornwall and R.E. Norton , ibid. 8, 3338 (1973). 

[45] D. H. Rischke, |nucl-th/0103050|. 

[46] M. Franz, Z. Tesanovic, |sond-niat/0012445 . 



12 



